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18 Abstract. A global climatology of tropospheric and stratospheric column ozone is 

19 derived by combining six years of Aura Ozone Monitoring Instrument (OMI) and 

20 Microwave Limb Sounder (MLS) ozone measurements for the period October 2004 

21 through December 2010. The OMI/MLS tropospheric ozone climatology exhibits large 

22 temporal and spatial variability which includes ozone accumulation zones in the tropical 

23 south Atlantic year-round and in the subtropical Mediterranean/Asia region in summer 

24 months. High levels of tropospheric ozone in the northern hemisphere also persist in 

25 mid-latitudes over the eastern North American and Asian continents extending eastward 

26 over the Pacific Ocean. For stratospheric ozone climatology from MLS, largest ozone 

27 abundance lies in the northern hemisphere in the latitude range 70°N-80°N in February- 

28 April and in the southern hemisphere around 40°S-50°S during months August-October. 

29 The largest stratospheric ozone abundances in the northern hemisphere lie over North 

30 America and eastern Asia extending eastward across the Pacific Ocean and in the 
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31 southern hemisphere south of Australia extending eastward across the dateline. With the 

32 advent of many newly developing 3D chemistry and transport models it is advantageous 

33 to have such a dataset for evaluating the performance of the models in relation to 

34 dynamical and photochemical processes controlling the ozone distributions in the 

35 troposphere and stratosphere. The OMI/MLS ozone gridded climatology data, both 

36 calculated mean values and RMS uncertainties are made available to the science 

37 community via the NASA total ozone mapping spectrometer (TOMS) website http://toms 

38 gsfc.nasa.gov. 

39 

40 1. Introduction. 

41 

42 In a previous paper Ziemke et al. [2006] combined ozone measurements from the Ozone 

43 Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS) onboard the Aura 

44 satellite to obtain global maps of tropospheric column ozone (TCO). The derivation of 

45 TCO was based upon a tropospheric ozone residual (TOR) method which involved 

46 subtracting MLS stratospheric column ozone (SCO) from OMI total column ozone after 

47 adjusting for calibration differences between the two instruments. The TOR concept, 

48 which was first applied by Fishman et al. [1990] involved measurements of total and 

49 stratospheric column ozone from two separate instruments on two separate satellites. 

50 Total column ozone was obtained from the Nimbus 7 TOMS UV backscatter instrument 

51 while SCO was obtained from Stratospheric Aerosols and Gas Experiment (SAGE) 

52 occultation instrument. Aside from calibration issues involving the use of two different 

53 satellite measurements, there was also a serious constraint in producing global data with 

54 adequate temporal and spatial coverage. Although TOMS total column ozone was daily 

55 with near global coverage, SAGE daily SCO measurements were limited to ~30 ozone 

56 profiles per day with about one month required to cover the latitude range 50°S-50°N. 

57 Chandra et al. [2003] later combined total column ozone from Nimbus 7 and Earth Probe 

58 TOMS with Upper Atmosphere Research Satellite (UARS) MLS ozone measurements to 

59 derive improved measurements of TCO extending from the tropics to the extra-tropics; 

60 that study was the first to use MLS ozone measurements to derive a long record of TCO 

61 spanning ~6 years. Having essentially coincident measurements along each orbital track 
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62 with current Aura OMI and MLS column ozone is a significant improvement from these 

63 previous studies in the ability to produce global maps of TCO from daily to longer 

64 timescales. 

65 

66 The study by Ziemke et al. [2006] used 2D interpolation of MLS to derive fields of SCO 

67 prior to deriving TCO. This method involved first applying a moving Gaussian 

68 interpolation to fill in missing MLS measurements along each daytime orbital track (on 

69 average about 14.6 south-to-north orbital tracks per day) followed by linear interpolation 

70 along longitude. This simple interpolation scheme for deriving daily gridded SCO fields 

71 works well provided that horizontal gradients in SCO are small. Meaningful daily maps 

72 of SCO can be determined from 2D interpolation excluding regions of polar night and the 

73 tropospheric wind jets. Other methods have been introduced for obtaining daily maps of 

74 SCO and TCO. Schoeberl et al. [2007] introduced a wind trajectory scheme using MLS 

75 ozone profiles and showed improvements for daily SCO and TCO in the extra-tropics. 

76 More recently Liu et al. [2010] describes an OMI-only ozone profile retrieval algorithm 

77 to derive daily fields of TCO and SCO. Other methods include data assimilation such as 

78 from the NASA Global Modeling and Assimilation Office (GMAO) [ Stajner , et al., 

79 2008] and also Aura TES ozone profile measurements [Zhang et al., 2010]. 

80 

81 All of these methods including simple 2D interpolation appear to do well in deriving 

82 global SCO and TCO for monthly averages, but accurate daily measurements remain 

83 generally constrained to excluding the wind jet regions. Even for regions where SCO 

84 does not vary appreciably, it remains to be shown that invoking more sophisticated 

85 methods beyond simple interpolation are an improvement because of potential induction 

86 of additional sources of error such as from assimilated wind fields or other ancillary data 

87 input to the algorithm. Tan et al. [2004] has shown that assimilated winds in the 

88 stratosphere have substantial errors caused by over-determination of mixing and 

89 entrainment rates in subtropical latitudes. Because of this, for a relatively long-lived 

90 constituent such as ozone in the lower stratosphere these errors propagate to errors in the 

91 tropical SCO field. It remains to be shown which of the various methods can provide 
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92 adequate accuracy and precision of SCO and TCO in the dynamical wind jet regions to 

93 derive true daily global maps (this is current collaborative work in progress). 

94 

95 The objective of this study is to develop a six- year global climatology of TCO and SCO 

96 from Aura OMI and MLS measurements. The importance of such a dataset is two-fold: 

97 (1) The dataset is useful as a benchmark test of basic seasonal-cycle and seasonal spatial 

98 variability present in photochemistry-transport models of tropospheric and stratospheric 

99 ozone, and also evaluating ozone precursors, and (2) the OMI/MLS dataset is also 

100 potentially useful as a priori information for ozone retrieval algorithms. It is noted that 

101 Martin et al. [2007] included an earlier version of OMI/MLS monthly mean TCO and 

102 derived an optimal assessment of global tropospheric nitrogen in the Goddard Earth 

103 Observing System Chemistry-transport model (GEOS-Chem) which was consistent with 

104 both the OMI/MLS ozone and the generated modeled ozone. In other words, by varying 

105 nitrogen in the model it was possible to yield consistent TCO between measurement and 

106 model, thus improving confidence in the assessed modeled nitrogen. 

107 

108 In the following, section 2 discusses the OMI and MLS satellite measurements, section 3 

109 compares the TCO with ozonesondes, sections 4 and 5 discusses TCO and SCO 

110 climatologies from OMI/MLS, and finally section 6 is a summary. 

111 

112 2. Aura OMI and MLS Ozone Measurements. 

113 

114 OMI and MLS are two out of a total of four instruments onboard the Aura spacecraft 

115 which is flown in a sun-synchronous polar orbit at 705 km altitude with a 98.2° 

116 inclination. Aura was launched in July 2004 and has been providing data measurements 

117 since August 2004 to the present. The spacecraft has an equatorial crossing time of 1:45 

118 pm (ascending node) with around 98.8 minutes per orbit (14.6 orbits per day on average). 

119 Schoeberl et al. [2006] provide an overview of the EOS Aura mission and discuss the 

120 various measurements from the four Aura instruments. 



122 OMI is a nadir-scanner which at visible (350-500 nm) and UV wavelength channels (UV- 

123 1: 270-314 nm; UV-2: 306-380 nm) detects backscattered solar radiance to measure 

124 column ozone with near global coverage (aside from polar night latitudes) over the Earth 

125 with a resolution of 13 km x 24 km at absolute nadir. Aside from ozone, OMI also 

126 measures Optical Centroid Cloud Pressure (OCCP), aerosols, NOa, SO 2 , HCHO, and 

127 several other trace gases in the troposphere and stratosphere [Levelt et al., 2006]. 

128 Measurements of ozone from OMI are determined using the OMTO3 v8.5 algorithm 

129 which is an extension of the TOMS v8 algorithm. A description of the TOMS v8 

130 algorithm may be obtained from the TOMS V8 CD DVD, or from the OMI Algorithm 

131 Theoretical Basis Document (ATBD) from the TOMS web page 

132 http://toms.gsfc.nasa.gov/version8/v8tonis_aibd.pdf). One main difference between the 

133 TOMS v8 and the OMT03 v8.5 algorithms is the treatment of clouds. The TOMS v8 

134 and earlier versions of 0MT03 use a cloud pressure climatology based on thermal 

135 infrared cloud-top pressures, whereas 0MT03 v8.5 uses OCCP derived with OMI by the 

136 rotational Raman scattering method. The use of simultaneously measured OCCP 

137 significantly improves estimates of total column ozone, especially in the presence of 

138 bright clouds [ Joiner and Vasilkov , 2006], The OMI instrument began having scan 

139 measurement errors beginning in year 2007 which are called row anomalies. The OMI 

140 instrument has 60 scan rows per orbit path and currently about 1/3 of the scan row 

141 measurements are affected. This problem is caused by partial blockage from material in 

142 the optical field of OMI but has not degraded appreciably since January 2009. 

143 

144 SCO is calculated for the OMI/MLS residual method using MLS v3.3 ozone profile 

145 measurements. The MLS instrument is a thermal-emission microwave limb sounder that 

146 measures vertical profiles of mesospheric, stratospheric, and upper tropospheric 

147 temperature, ozone, and several other important constituents such as CO and H 2 O from 

148 limb scans taken in the direction ahead of the Aura satellite orbital track. The MLS 

149 profile measurements are made about 7 minutes before OMI views the same location 

150 during ascending (daytime) orbital tracks. These we refer to as "co-located" 

151 measurements between OMI and MLS. MLS also measures ozone and other atmospheric 

152 constituents for descending nighttime orbits which on a given day can be up to +12 hours 
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153 different in time from OMI daytime measurements. With combined ascending and 

154 descending nodes MLS makes around 3500 vertical profile measurements over the Earth 

155 per day. This study includes only the ascending orbit co-located data from MLS for 

156 deriving SCO. Details regarding the instrument including spectrometers, spectral 

157 channels, calibration, and other topics are discussed by Waters et al. [2006J and in related 

158 papers in the same journal. Froidevaux et al. [ 2008 1 provides validation results for MLS 

159 v2.2 measurements of ozone and other trace gases. At the present time both v2.2 and 

160 v3.3 are provided to the science community. While v2.2 retrieval has 37 pressure levels, 

161 v3.3 has 55 pressure levels and other improvements; however, v3.3 also has more outliers 

162 and missing data in the ozone profile measurements than with v2.2. Our analysis of SCO 

163 from MLS show's that there is little difference between using v2.2 or v3.3 other than a 

164 small systematic offset (v3.3 minus v2.2) at all latitudes of about +2.5 DU. As was 

165 similarly done by Ziemke et al. [2006], MLS v3.3 SCO was adjusted to CCD SCO from 

166 OMI by subtracting this offset from MLS SCO. Information regarding the MLS v3.3 

167 ozone measurements including a data quality description document is available online 

168 from the NASA Data and Information Services Center webpage 

169 (http://disc.sci.gsfc.nasa.gov/gesNews/mls_new_data_version_release). 

170 

171 For the OMI/MLS residual method SCO is determined from vertically integrated MLS 

172 ozone profiles which is subtracted from OMI total column ozone to derive TCO. 

173 Tropopause pressure, which separates tropospheric from stratospheric column ozone 

174 comes from National Centers for Environmental Prediction (NCEP) using the World 

175 Meteorological Organization (WMO) 2K-km 1 lapse rate tropopause definition. SCO 

176 from MLS is determined by pressure integration of ozone volume mixing ratio profiles 

177 from 0.0215 hPa down to the NCEP tropopause. The MLS ozone profile measurements 

178 were linearly interpolated in log-pressure to the NCEP tropopause pressure to derive 

179 SCO. MLS SCO (in Dobson Units, DU; 1 DU = 2.69xl0 20 molecules-m' 2 ) was 

180 determined by standard log-pressure integration of ozone partial pressure: 

rPtropopause 

181 SCO = 0.79 1 XP dlnP , where X is ozone volume mixing ratio in units ppbv 

+).0215 hPa 

182 and P is pressure in units hPa. The recommended range for scientific analysis of MLS 
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183 v3.3 ozone profiles is 0.0215-261 hPa. As was done by Ziemke et al. [2006], nearly 

184 global SCO from MLS for each day was achieved by including ozone retrievals down to 

185 316 hPa. The uncertainty in derived SCO from MLS by including the 316 hPa level 

186 beyond 261 hPa can be estimated. Documentation for MLS v3.3 ozone measurements 

187 suggests RMS uncertainties of about 0.03 ppmv at 261 hPa and 0.05 ppmv at 316 hPa. 

188 Using the above integration formula, an upper-bound estimate of the RMS uncertainty in 

189 261-316 hPa column ozone from MLS is then about 3.5 DU. 

190 MLS SCO data were initially binned to 1° latitude x 1.25° longitude to be compatible 

191 with OMI level-3 (L3) gridded total column ozone. Tropopause pressures from NCEP 

192 analyses were re-binned to this same resolution from a coarser 2.5° x 2.5° gridding. It is 

193 noted for MLS limb measurements that the horizontal optical path is about 300 km which 

194 is larger than the horizontal size of OMI L3 gridded data, but is comparable to the size of 

195 original NCEP gridded measurements. To derive a high density SCO field we have used 

196 the two-step spatial interpolation of Ziemke et al. [2006]. The interpolation for SCO 

197 includes first a moving 2D (latitude/longitude) Gaussian window along daytime orbit to 

198 fill in intermittent gaps along-track for MLS SCO, followed secondly by a linear 

199 interpolation along longitude between MLS SCO data. This interpolation approach 

200 preserves the along-track measurements of SCO from MLS at all latitudes. NCEP 

201 measurements of tropopause pressure were re-binned to the same 1° latitude x 1.25° 

202 longitude resolution. Following derivation of daily maps of SCO and TCO at l°x 1.25° 

203 resolution, the data were averaged monthly in 5° x 5° and 10° x 10° bins. We have 

204 evaluated OMI and MLS ozone data for the time period I October 2004 - 3 1 December 

205 2010 which represents approximately six years of continuous measurements. The TCO 

206 and SCO monthly climatology fields were smoothed along longitude using Fourier 

207 analysis by retaining zonal wave numbers 0-12 which both reduces noise and ensures 

208 periodic continuation across the dateline. 

209 

210 3. Comparisons of QMI/MLS Tropospheric Ozone with Ozonesondes. 

211 


7 



212 Validation of OMI/MLS TCO is based upon Southern Hemisphere Additional 

213 OZonseondes (SHADOZ) and World Ozone and Ultraviolet radiation Data Center 

214 (WOUDC) ozonesonde measurements whereby monthly means from OMI/MLS are 

215 compared with monthly ensemble averages from the ozonesondes. The ozonesonde data 

216 represent for each station an average of all existing ozonesonde measurements in a given 

217 month (which could vary from one to several). OMI/MLS TCO was converted to ozone 

218 mean volume mixing ratio and then compared with mean volume mixing ratio for both 

219 the SHADOZ and WOUDC ozonesondes. Ozone mean volume mixing ratio for the 

220 tropical SHADOZ stations involved pressure averaging of ozone for ground-to-120 hPa 

221 while for extra-tropical WOUDC stations this was ground-to-350 hPa. The ozonesonde 

222 data for our comparisons extend for years 2004-2009 for SHADOZ and 2005-2008 for 

223 WOUDC. Of the several SHADOZ and WOUDC stations processed, 39 of them had 

224 measurements that overlapped with the OMI/MLS time period. 

225 

226 Figure 1 and Figure 2 show tropospheric ozone time series comparisons for several 

227 stations from SHADOZ and WOUDC, respectively. Also included in these figures are 1 - 

228 1 scatter diagrams of the combined data. The stations in Figures 1 and 2 were chosen not 

229 just because these are common stations but also because of large annual-cycle variability 

230 (20-30 ppbv peak-to-peak) with nearly complete month-to-month temporal coverage. 

231 

232 There are important issues involved when comparing ozonesondes with satellite ozone 

233 measurements. Balloon ozonesondes represent wind advected measurements along 

234 trajectory paths and are very different than the satellite measurements which provide 

235 ozone averages over broad regions. Also, ozone in the troposphere (such as shown by 

236 Avery et al. [2010]) exhibits large changes of 30-50 ppbv over horizontal distances of 

237 only a few km or less. With such enormous spatial variability the most one can gain in 

238 comparing satellite measured ozone versus ozonesondes is to evaluate seasonal 

239 variability and climatological means, and extreme cases such as inter-annual variability 

240 events in the tropical Pacific caused by ENSO. Table 1 provides a list of comparisons 

241 between OMI/MLS and the SHADOZ/WOUDC ozonesondes with station latitudes 

242 arranged from northern-most at top to southern-most at bottom. Listed for each station is 



243 the total number of daily profile measurements, mean values for OMI/MLS and 

244 ozonesondes, and standard RMS of their differences. For most of the stations listed in 

245 Table l the RMS values lie between 4 ppbv and 10 ppbv. For stations northward of 50°N, 

246 RMS is greater than 10 ppbv except for Debilt and Valentia. The larger RMS numbers in 

247 NH high latitudes are partly explained by the sparse nature of monthly measurements 

248 from both OMI/MLS and ozonesondes. Figure 3 compares climatological means for 

249 OMI/MLS and ozonesondes. The left panel in Figure 3 represents station latitudes 25°S- 

250 50°N and the right panel is the same but includes stations pole ward of 50°N. 

251 

252 4. Comparisons of OMI/MLS and LLM Tropospheric Ozone Climatology. 

253 

254 Logan [1999] provided an extensive analysis of global ozone based upon ozonesonde 

255 measurements. More recently McPeters et al. [2007] expanded the ozonesonde 

256 evaluation of Logan [1999] and derived a global zonal mean climatology of total 

257 atmospheric ozone as a function of latitude, altitude, and month of year. This 

258 climatology was determined by combining ozonesondes with satellite ozone 

259 measurements from SAGE and MLS. The final climatology product is referred to as the 

260 Labow-Logan-McPeters (LLM) climatology and is currently used in the OMT03 v8.5 

261 algorithm processing for both OMI and TOMS ozone retrievals. 

262 

263 Figure 4 shows a comparison of LLM TCO climatology (top panel) with OMI/MLS TCO 

264 climatology (middle panel) and their difference (bottom panel). Fundamental 

265 characteristics are generally consistent between the datasets regarding seasonal cycles 

266 including the transition in peak TCO in the NH from tropics to mid-latitudes when going 

267 from spring to summer months. LLM minus OMI/MLS TCO difference is larger by 

268 several DU to more than +5 DU for some months and latitudes in the subtropics and mid- 

269 latitudes. These differences are largely explained as an ozonesonde station location issue 

270 as the LLM climatology is determined from TCO measurements over land whereas TCO 

271 from OMI/MLS in Figure 4 is averaged along all longitudes. 

272 
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273 That is, because TCO over ocean is generally smaller than over land (as discussed later 

274 for Figure 5), OMI/MLS zonal mean TCO will be smaller than LLM TCO. In contrast, 

275 in the high latitudes of both hemispheres in Figure 4 the offset differences are instead 

276 negative exceeding -5 DU. These differences in the higher latitudes may be explained by 

277 both a sparse number of ozonesonde measurements and also high solar-zenith angles 

278 which can produce substantial errors in the OMI and/or MLS ozone retrievals. 

279 

280 5. OMI/MLS Global TCO Climatology Maps. 

281 

282 The first study to derive a global TCO climatology maps from satellite measurements 

283 was Fishman et al. [1990J. As noted in the Introduction, tropospheric ozone determined 

284 from TOMS/SAGE residual was limited both temporally and spatially because of sparse 

285 SAGE ozone profile measurements. An improvement was made for the SAGE sampling 

286 problem by Fishman et al. [2003] by using instead NOAA SBUV to derive SCO, 

287 however there is little ozone profile information obtained from SBUV below the ozone 

288 number density peak (i.e., atmospheric pressures greater than -30-40 hPa). To alleviate 

289 this SBUV ozone profile problem Fishman et al. [2003] applied an empirical correction 

290 to SBUV ozone profiles to improve the SCO fields. It is anticipated that the use of Aura 

291 MLS ozone profile measurements which extend down to 261 hPa for v3.3 substantially 

292 improves global climatology estimates of SCO and also TCO. 

293 

294 A six-year TCO climatology from OMI and MLS ozone is shown in Figure 5(a-c) for 

295 each of the months January-December. Figure 5 for illustration plots the 10°xl0 H dataset. 

296 TCO in high latitudes is flagged as missing, both for obvious reason involving polar night 

297 latitudes (i.e., where there are no OMI measurements) and additional subjective data 

298 flagging. Subjective data flagging included OMI scenes with solar zenith angles greater 

299 than 82° or questionable derived TCO product values in high latitudes. 

300 

301 The main features of the TCO climatology in Figure 5 are briefly summarized. In 

302 tropical latitudes TCO for each month is characterized by low amounts in the Pacific 

303 -15-20 DU with much higher columns -35-45 DU in the Atlantic. Lowest TCO in the 



304 tropics lies in the western Pacific in July-September with numbers less than 15 DU. 

305 These low columns are largely a manifestation of deep convection and vertical injection 

306 of low marine boundary layer/low tropospheric ozone into the middle and upper 

307 troposphere. The highest TCO in the tropics occurs in September-October in the 

308 southern Atlantic region extending eastward toward Australia with numbers 40-45 DU. 

309 These annually recurring high values in the southern hemisphere are associated with 

310 planetary-scale transport of ozone rich air mass and ozone sources including effects from 

311 STE and biomass burning. TCO in the northern hemisphere has large topographic 

312 variability year-round, with largest columns occurring around May-July over the eastern 

313 US eastward over the Atlantic Ocean toward Europe, Mediterranean/westem Asia, and 

314 eastern Asia extending eastward over the Pacific Ocean toward North America. 

315 

316 Figure 6 shows line plots of TCO as a function of longitude for 30°-40°N and 0°-10°S. 

317 These latitude bands and months were chosen to show several large annually recurring 

318 features present in global TCO. The line plots in Figure 6 also include ±2o uncertainties 

319 for illustration where a is calculated standard RMS error of the mean. TCO in these two 

320 latitude bands has large spatial variability. The left panel shows zonal variability caused 

321 by a Mediterranean- Asian accumulation region (centered about 30°E) and the mountains 

322 of the Tibetan Plateau (centered about 90°E). The right panel shows large contrast in 

323 TCO from low values in the tropical Pacific and high values in the tropical Atlantic (i.e., 

324 about 20 DU versus 40 DU, respectively). The 5°x5° and 10°xl0° monthly TCO 

325 climatology fields along with their RMS uncertainty fields can be obtained via link from 

326 the TOMS webpage (http://torns.gs fc . nasa . gov) . 

327 

328 For applications such as evaluating 3D chemistry-transport models of the atmosphere it is 

329 useful to derive a zonal mean climatology of TCO from OMI/MLS. A simple model 

330 comparison with zonal mean TCO climatology can reveal fundamental offsets or annual- 

331 cycle differences that can be used to aid in the development of models. Ziemke et al. 

332 [2006] compared an earlier version of zonal mean TCO between OMI/MLS and the 

333 Global Modeling Initiative (GMI) 3D model and showed good agreement in latitude 

334 gradients and temporal variability except for a 5-10 DU offset in the NH extending from 



335 the subtropics out to mid-latitudes. It was surmised that errors in OMI from aerosols 

336 could be a contributing factor for the offset, but the primarily indication was that the 

337 model was over-determining TCO by several DU in the NH subtropics/mid- latitudes in 

338 most months. Figure 7 shows the six-year climatology of zonal mean TCO. Lowest 

339 TCO occurs in the SH tropics around 1 0°S around January- April (<24 DU) and also in 

340 the NH and SH mid-latitudes around late fall to winter-spring (<30 DU). Largest TCO 

341 occurs in the NH mid-latitudes in June-July (>42 DU) and in the SH subtropics during 

342 September-November (>39 DU). There is a transition in the NH from peak values in 

343 spring months (March-May) in the tropics/subtropics to peak values in summer (June- 

344 July) in the mid-latitudes. This pattern shift from spring to summer with latitude is 

345 caused largely by a coupling of spring-summer STE with ozone produced from pollution 

346 events in summer months. It is noted that the GMI model evaluated by Ziemke et al. 

347 [2006J showed all the basic temporal and spatial features in Figure 7 despite an offset 

348 difference in the NH mid- latitudes. 

349 

350 6. MLS Global SCO Climatology Maps. 

351 

352 An SCO climatology from MLS is shown in Figure 8(a-c) for each of the months 

353 January-December. As in Figure 4 the data binning is at 10°xl0° resolution which for 

354 SCO provides nearly global coverage out to central latitudes ±85°. The main features 

355 include a contrast between small zonal variability of SCO in the tropics and large 

356 planetary-scale variability in the middle and high latitudes in both hemispheres. SCO in 

357 Figure 8 should be interpreted carefully in the middle and high latitudes particular to 

358 these years 2004-2010 since SCO patterns in both hemispheres can exhibit substantial 

359 inter-annual differences in zonal variability caused by stratospheric sudden warming 

360 events and the morphology of the breakup of the middle atmosphere polar vortex in the 

361 NH and SH. Zonal mean SCO from MLS is shown in Figure 9. SCO in Figure 9 is 

362 largest in spring months in both hemispheres. In the Northern hemisphere SCO is largest 

363 around February-March at high latitudes and in the Northern Hemisphere SCO is largest 

364 around September-October in mid-latitudes. The maximum SCO in the NH occurs in 

365 conjunction with high ozone over the central Asian continent extending eastward across 
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366 North America as shown in Figure 7a. SCO in the SH maximizes equator-ward of the 

367 polar vortex in September-October and as Figure 7c shows, originates primarily from the 

368 Pacific about the dateline. Table 3 shows the calculated zonal mean SCO climatology 

369 values given in 5° latitude bands. 

370 

371 7. Summary. 

372 

373 A six-year global climatology of tropospheric column ozone (TCO) and stratospheric 

374 column ozone (SCO) is derived from combining Aura OMI and MLS measurements for 

375 the period October 2004-December 2010. The climatology is useful for several purposes 

376 including evaluation of 3D chemistry-transport models and self-generating global 

377 circulation models of the atmosphere to assess model performance. By comparing basic 

378 seasonal cycles and seasonally-varying spatial variability in TCO and SCO in models 

379 with the climatology one can identify and correct problems in models regarding errors in 

380 wind fields or invoked photochemical reactions/rates, and errors in emissions including 

381 ozone precursors. Another useful application for the global climatology of TCO and 

382 SCO includes producing a climatology of radiative forcing coming from upper and lower 

383 atmospheric ozone. The total column ozone climatology determined from adding 

384 together SCO and clear-sky TCO can also generate a global clear-sky surface UV 

385 climatology as a function of region and time of year. 

386 

387 Algorithms for retrieving ozone such as from remote sensing satellite instruments rely on 

388 assumed a priori information of the ozone concentrations. The current 0MT03 

389 algorithm invokes the Labow-Logan-McPeters (LLM) zonal mean ozone climatology 

390 [McPeters, et al„ 2007J which is determined by combining ozonesondes with SAGE and 

391 Aura MLS ozone profile measurements. The LLM climatology is given as a function of 

392 month, latitude, and altitude. A future task is to combine the LLM and the OMI/MLS 

393 climatology products to generate a single 3D (latitude, longitude, altitude) 12-month 

394 climatology. 

395 
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396 The TCO and SCO climatology data from OMI/MLS are provided to the community in 

397 ASCII formatted tables with IDL and Fortran readers as a link from the TOMS homepage 

398 (http://toms.gsfc.nasa.gov). The measurements along with their computed RMS 

399 uncertainties are given in these tables at 5° x 5° and 10° x 10° (latitude x longitude) 

400 resolution. 

401 
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477 Table and Figure captions 

478 

479 Table 1. Statistical comparisons between OMI/MLS and SHADOZ/WOUDC sonde 

480 tropospheric mean volume mixing ratio (in ppbv) with the station latitudes arranged from 

481 northern-most at top to southern-most at bottom. Included for each station is the total 

482 number of daily profde measurements, mean values for OMI/MLS and ozonesondes, and 

483 standard RMS of their differences. 

484 

485 Table 2. Tropospheric column ozone zonal mean climatology (in Dobson Units) from 

486 OMI/MLS at 5° latitude resolution. 
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487 

488 Table 3. Global stratospheric column ozone zonal mean climatology (in Dobson Units) 

489 derived from MLS integrated ozone profiles at 5° latitude resolution. 

490 

491 Figure 1. Comparison time series of tropospheric ozone mean volume mixing ratio (in 

492 ppbv) between OMI/MLS (dark curves) and several selected common station 

493 measurements from SHADOZ (light curves). Also included at far right is a 1-1 

494 comparison scatter plot for the measurements. 

495 

496 Figure 2. Similar to Figure 1 but for selected extra-tropical WOUDC ozonesonde 

497 stations. 

498 

499 Figure 3. (left) Scatter plots of ozonesonde versus OMI/MLS mean tropospheric ozone 

500 for ozonesonde station locations lying between 25°S and 50°N. (right) Same as left panel, 

501 except for extended latitude range 25°S to 90°N. The values plotted are time series 

502 averages in units ppbv. 

503 

504 Figure 4. (top) Labow-Logan_McPeters (LLM) tropospheric column ozone climatology. 

505 (middle) OMI/MLS tropospheric column ozone climatology, (bottom) LLM minus 

506 OMI/MLS climatology difference. All measurements are in Dobson Units. The colors in 

507 the panels going from blue/black to red represent smallest to largest (or most positive) 

508 values, respectively. 

509 

510 Figure 5a. Tropospheric column ozone climatology (in Dobson Units) for months 

511 January-April from OMI/MLS residual ozone measurements. The colors in the panels 

512 going from blue to red represent smallest to largest values, respectively. 

513 

514 Figure 5b. Same as Figure 5a but for months May-August. 

515 

516 Figure 5c. Same as Figure 5b but for months September-December. 

517 
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Figure 6. Line plots of tropospheric column ozone (in Dobson Units) as a function of 
longitude for the latitude bands 30°-40°N and 0°-10°S. The plots include ±2o 
uncertainties for illustration where a is the calculated standard RMS error of the mean. 

Figure 7. Zonal mean climatology of tropospheric column ozone (in Dobson Units) 
derived from October 2004 - December 2010 OMI/MLS at 5° latitude resolution. The 
shading going from dark to light represents smallest to largest values, respectively. 

Figure 8a. Similar to Figure 5a but for stratospheric column ozone (in Dobson Units). 
The colors in the panels going from blue to red represent smallest to largest values, 
respectively. 

Figure 8b. Similar to Figure 5b but for stratospheric column ozone (in Dobson Units). 

Figure 8c. Similar to Figure 5c but for stratospheric column ozone (in Dobson Units). 

Figure 9. Zonal mean climatology of stratospheric column ozone (in Dobson Units) 
from October 2004 - December 2010 MLS at 5° latitude resolution. The shading going 
from dark to light represents smallest to largest values, respectively. 


Tables 


Table 1. Statistical comparisons between OMI/MLS and SHADOZ/WOUDC sonde 
tropospheric mean volume mixing ratio (in ppbv) with the station latitudes arranged from 
northern-most at top to southern-most at bottom. Included for each station is the total 
number of daily profde measurements, mean values for OMI/MLS and ozonesondes, and 
standard RMS of their differences. 


Station 

N 

OMI/MLS 

Sonde 

Diff 

RMS 

Alert(82.5N,62.3W) 

27 

60.2 

49.4 

19.4 

Eureka(80.QN,85.9W) 

22 

56.9 

50.5 

13.3 
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Resolute(74.7N,95.0W) 

33 

61.9 

49.7 

19.1 

Neumayer(70.7S,8.3W) 

26 

38.4 

30.8 

17.7 

Lerwick(60. 1 N, 1 .2W) 

46 

57.2 

55.2 

10.7 

Churchill(58.7N,94. 1 W) 

36 

61.1 

50.6 

14.3 

StonypIain(53.6N,l 14. 1W) 

45 

53.1 

48.2 

11.0 

Goosebay(53.3N,60.4W) 


57.9 

51.4 

13.2 

Legionwo(52.4N,2 1 .OE) 


51.8 

56.5 

12.1 

Lindenberg(52.2N, 14. IE) 


51.0 

55.9 

11.1 

Debilt(52. 1N,5.2E) 


52.1 

54.7 

8.4 

Valentia(5 1 ,9N, 10.3W) 





Brattslake(50.2N, 104.7W) 





Praha(50.()N, 14.5E) 


50.4 

54.1 

7.9 

Kelowna(49.9N,l 19.4W) 

36 

51.1 

51.8 

12.4 

Payerne(46.8N,7.0E) 

48 


52.2 

6.6 

Egbert(44.2N,79.8W) 

41 

58.3 

57.8 

9.3 

Barajas(40.5N,3.6W) 

47 

55.9 

53.9 

7.0 

Boulder(40.3N,105.2W) 

36 

54.0 

56.0 

6.6 

Ankara(40.0N,32.9E) 

43 

57.5 

59.2 

9.9 

Wallops(37.9N,75.5W) 





Huntsville(34.7N,86.6W) 




12.4 

Isfahan(32.5N,5 1 ,7E) 




11.5 

Hongkong(22.3N,l 14.2E) 


49.7 

50.9 

7.8 

Hilo( 1 9.4N, 1 55. OW) 



50.9 


Alajuela( 1 0.0N,84.2W) 


■ 

40.5 

4.2 

Heredia) 10.0N, 84. 1W) 


ym 

43.5 

5.9 

Panama)7.8N,80.3W) 

m 

38.7 

44.0 

6.5 

Kuala(2.7N, 1 0 1 ,7E) 

El 

34.2 

35.2 

5.6 

Sancr(0.9S,89.6W) 


38.4 

36.7 

5.0 

Nairobi) l.3S,36.8E) 

52 

37.8 

46.2 

9.9 

Malindi(3.0S,40,2E) 

13 

45.0 

51.3 

7.6 
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Natal(5.4S,35.4W) 

58 

47.7 

49.0 

6.0 

Java(7.5S,l 12.6E) 

44 

35.3 

34.8 

6.7 

Ascen(8.0S, 1 4.4W) 

56 

50.9 

55.2 

6.7 

Samoa) 14.2S, 170. 6W) 

56 

31.5 

30.8 

4.5 

Fiji) 1 8. 1 S, 1 78.4E) 

22 

34.7 

34.2 

4.7 

Reunion(2 1 .0S.55.5E) 

52 

50.5 

59.9 

11.7 


547 

548 

549 Table 2. Tropospheric column ozone zonal mean climatology (in Dobson Units) from 

550 OMI/MLS at 5° latitude resolution. 

551 


Latitudes 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

60N-65N 

35.8 

34.6 

33.1 

33.7 

34.3 

33.8 

33.8 

30.7 

28.2 

27.9 

33.8 

35.8 

55N-60N 

33.0 

32.2 

33.0 

33.7 

34.6 

35.6 

35.5 

32.9 

29.5 

28.1 

30.3 

33.6 

50N-55N 

31.1 

32.1 

32.9 

33.5 

35.2 

36.7 

36.9 

34.7 

30.8 

28.4 

28.8 

30.1 

45N-50N 

29.4 

30.8 

31.6 

33.7 

35.9 

38.2 

39.1 

37.4 

33.2 

28.9 

27.9 

28.0 

40N-45N 

27.6 

28.9 

30.9 

34.1 

37.0 

40.9 

42.2 

40.7 

35.7 

30.2 

27.3 

27.3 

35N-40N 

27.0 

27.7 

30.9 

35.2 

38.9 

43.0 

44.0 

42.3 

37.9 

32.4 

28.4 

27.3 

30N-35N 

28.2 

28.6 

32.9 

38.0 

41.2 

42.9 

42.3 

40.4 

37.3 

33.5 

30.9 

29.0 

25N-30N 

30.2 

31.6 

35.9 

40.5 

42.0 

39.8 

38.3 

36.8 

35.0 

33.1 

31.9 

30.6 

20N-25N 

31.6 

32.4 

36.9 

40.5 

40.4 

37.1 

35.3 

33.6 

33.0 

32.0 

31.7 

31.6 

15N-20N 

31.1 

31.3 

35.2 

38.0 

37.1 

34.3 

32.0 

30.0 

29.8 

29.7 

30.2 

31.0 

10N-15N 

29.2 

29.3 

32.7 

34.1 

32.6 

30.3 

28.2 

26.8 

27.3 

27.7 

28.6 

29.5 

5N-10N 

27.2 

27.6 

29.9 

29.8 

27.7 

26.5 

25.8 

25.1 

26.1 

26.3 

26.6 

27.7 

0-5N 

25.2 

25.5 

26.8 

26.0 

24.9 

25.3 

25.8 

26.4 

27.2 

26.9 

26.4 

26.4 

0-5S 

24.4 

24.0 

24.8 

24.1 

25.0 

26.9 

27.5 

28.7 

30.1 

29.9 

28.5 

26.9 

5S-10S 

24.2 

22.9 

23.4 

23.5 

25.4 

27.8 

28.5 

30.0 

32.0 

32.3 

30.4 

27.7 

10S-15S 

24.5 

23.0 

23.6 

23.9 

26.0 

28.7 

29.5 

30.8 

33.3 

33.9 

32.1 

28.7 

15S-20S 

26.6 

24.7 

24.8 

25.2 

27.1 

29.8 

31.0 

32.5 

35.5 

36.4 

34.2 

30.5 

20S-25S 

29.8 

27.9 

27.6 

27.6 

28.7 

30.9 

32.7 

35.1 

38.2 

39.5 

37.3 

33.5 

25S-30S 

32.4 

30.6 

29.9 

29.3 

29.8 

31.4 

34.2 

36.8 

40.2 

41.2 

39.3 

36.2 


20 





30S-35S 

33.9 

31.8 

30.5 

28.6 

28.3 

30.0 

33.4 

36.2 

38.9 

39.5 

37.5 

35.9 

35S-40S 

32.2 

30.4 

28.4 

26.4 

26.2 

27.6 

30.7 

33.9 

35.4 

35.0 

32.5 

31.7 

40S-45S 

28.0 

27.2 

26.0 

25.6 

25.6 

26.5 

28.5 

30.5 

31.6 

30.3 

28.3 

27.1 

45S-50S 

25.6 

25.8 

25.7 

25.9 

25.8 

25.9 

26.7 

28.2 

30.6 

29.9 

28.0 

25.8 

50S-55S 

25.3 

26.5 

27.8 

27.5 

27.8 

29.7 

28.5 

28.0 

31.4 

30.7 

28.4 

25.4 

55S-60S 

23.6 

26.5 

28.4 

26.8 

28.6 

28.5 

29.4 

28.8 

30.2 

29.5 

26.7 

23.7 

60S-65S 

22.6 

25.8 

27.5 

26.6 

29.8 

33.2 

33.6 

28.3 

27.0 

31.2 

25.8 

22.0 


552 

553 

554 Table 3. Global stratospheric column ozone zonal mean climatology (in Dobson Units) 

555 derived from MLS integrated ozone profdes at 5° latitude resolution. 

556 


Latitudes 

Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

85N-90N 

311 

350 

352 

340 

324 

297 

269 

253 

243 

246 

267 

279 

80N-85N 

314 

363 

361 

359 

330 

300 

269 

252 

244 

251 

270 

281 

75N-80N 

319 

363 

363 

359 

330 

301 

272 

254 

248 

254 

275 

287 

70N-75N 

326 

359 

361 

359 

333 

303 

275 

262 

256 

259 

279 

293 

65N-70N 

330 

359 

359 

355 

332 

304 

280 

269 

262 

265 

280 

299 

60N-65N 

333 

359 

359 

352 

331 

306 

286 

274 

267 

269 

283 

306 

55N-60N 

332 

355 

354 

347 

329 

307 

291 

277 

270 

270 

282 

309 

50N-55N 

329 

349 

348 

340 

325 

306 

290 

276 

267 

267 

279 

308 

45N-50N 

321 

338 

338 

331 

318 

299 

280 

269 

262 

260 

273 

302 

40N-45N 

307 

320 

324 

318 

304 

286 

266 

259 

254 

252 

265 

289 

35N-40N 

284 

292 

300 

298 

286 

271 

256 

252 

248 

243 

252 

269 

30N-35N 

257 

262 

272 

275 

269 

259 

251 

248 

244 

237 

237 

247 

25N-30N 

234 

238 

248 

255 

256 

252 

248 

246 

242 

234 

229 

230 

20N-25N 

221 

224 

234 

242 

246 

246 

245 

244 

240 

232 

224 

219 

15N-20N 

215 

218 

227 

235 

240 

242 

244 

245 

241 

232 

223 

215 

10N-15N 

214 

216 

224 

232 

237 

240 

244 

245 

242 

232 

223 

214 

5N-10N 

215 

217 

223 

230 

234 

237 

240 

242 

240 

231 

224 

216 

0-5N 

219 

221 

226 

231 

232 

233 

236 

238 

237 

229 

225 

219 
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0-5S 

223 

224 

228 231 

230 

229 

5S-10S 

225 

226 

228 229 

227 

224 

10S-15S 

227 

227 

227 227 

223 

221 

15S-20S 

228 

226 

226 225 

222 

222 

20S-25S 

228 

225 

225 225 

224 

226 

25S-30S 

229 

226 


234 

30S-35S 

232 

230 


249 


35S-40S 


40S-45S 


229 230 233 233 228 225 222 

224 225 228 230 228 227 225 ~ 

221 223 226 230 229~1>30 228~ 

222 225 228 234 234 233 230~ 

226 230 235 242 241 239 232~ 


60S-65S I 274 I 264 
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QMI/MLS and Ozonesondes in the Tropics 


591 

592 

593 

594 

595 

596 

597 


Ascert Lat:~7.98 Lon:- 14.42 



Month 


Natal Lat:-5.42 Lon:-35.38 



Month 


Hilo Lat: 19.40 Lon:~155.Q 



J04 JOS J06 J07 J08 JOS JlO 

Month 


Figure 1. Comparison time series of tropospheric ozone mean volume mixing ratio (in 
ppbv) between OMI/MLS (dark curves) and several selected common station 
measurements from SHADOZ (light curves). Also included at far right is a 1-1 
comparison scatter plot for the measurements. 


OMI/MLS and Ozonesondes in the Extra-Tropics 


598 

599 

600 

601 

602 


Ankara Lat;39.97 Lon;32.86 





J06 J07 JOS JOS 

Month 

Wallops Lat:37.93 Lon:-75.47 


J07 

Month 



OMI/MLS and Sonde Tropo 03 



Figure 2. Similar to Figure 1 but for selected extra-tropical WOUDC ozonesonde 
stations. 
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Tropo 03 Mean Values 25S-50N 



Tropo 03 Mean Values 25S-9QN 



Figure 3. (left) Scatter plots of ozonesonde versus OMI/MLS mean tropospheric ozone 
for ozonesonde station locations lying between 25°S and 50°N (right) Same as left panel, 
except for extended latitude range 25°S to 90°N. The values plotted are time series 
averages in units ppbv. 
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612 Figure 4. (top) Labow-Logan_McPeters (LLM) tropospheric column ozone climatology. 

613 (middle) OMI/MLS tropospheric column ozone climatology, (bottom) LLM minus 

614 OMI/MLS climatology difference. All measurements are in Dobson Units. The colors in 

615 the panels going from blue/black to red represent smallest to largest (or most positive) 

616 values, respectively. 

617 

618 
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Tropospheric 03 (DU) January 
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180 -120 -60 0 60 120 
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Figure 5a. Tropospheric column ozone climatology (in Dobson Units) for the months 
January-April from OMI/MLS residual ozone measurements. The colors in the panels 
going from blue to red represent smallest to largest values, respectively. 
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Tropospheric 03 (DU) June 
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Tropospheric 03 (DU) July 
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Tropospheric 03 (DU) August 


Longitude 

Figure 5b. Same as Figure 5a but for months May 
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Tropospheric 03 (DU) September 



180 -120 -60 0 60 120 180 
Longitude 
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Tropospheric 03 (DU) December 
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Figure 5c. Same as Figure 5b but for months September-December. 
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Figure 6. Line plots of tropospheric column ozone (in Dobson Units) as a function of 
longitude for the latitude bands 30°-40°N and 0°-10°S. The plots include ±2o 
uncertainties for illustration where a is the calculated standard RMS error of the mean. 
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OMI/MLS Tropo Ozone (DU) Zonal Means 


Figure 7. Zonal mean climatology of tropospheric column ozone (in Dobson Units) 
derived from October 2004 - December 2010 OMI/MLS at 5° latitude resolution. The 
shading going from dark to light represents smallest to largest values, respectively. 


Stratospheric 03 (DU) January 


Stratospheric 03 (DU) February 


Longitude 

Stratospheric 03 (DU) March 
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Stratospheric 03 (DU) April 


Longitude Longitude 

Figure 8a. Similar to Figure 5a but for stratospheric column ozone (in Dobson Units). 
The colors in the panels going from blue to red represent smallest to largest values, 
respectively. 
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Figure 8b. Similar to Figure 5b but for stratospheric column ozone (in Dobson Units) 









Stratospheric 03 (DU) September 


Stratospheric 03 (DU) October 


Longitude 

Stratospheric 03 (DU) November 
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Stratospheric 03 (DU) December 
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Longitude 


Figure 8c. Similar to Figure 5c but for stratospheric column ozone (in Dobson Units) 
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656 Figure 9. Zonal mean climatology of stratospheric column ozone (in Dobson Units) 

657 from October 2004 - December 2010 MLS at 5° latitude resolution. The shading going 

658 from dark to light represents smallest to largest values, respectively. 
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